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ABSTRACT 

We combine GALEX (ultra-violet; UV) and SDSS (optical) photometry to study the 
recent star formation histories of ^ 100 field galaxies on the optical red sequence, a 
large fraction of which exhibit widespread signs of disturbed morphologies in deep 
optical imaging that are consistent with recent merging events. More than 70% of 
bulge-dominated galaxies in this sample show tidal features at a surface brightness 
limit of /^ ^ 28 mag arcsec"^. We find that, while they inhabit the optical red sequence, 
they show a wide spread in their UV colours (~ 4 mags), akin to what has been 
discovered recently in the general early-type population. A strong correlation is found 
between UV colour and the strength of the tidal distortions, such that the bluest 
galaxies are more distorted. This strongly suggests that the blue UV colours seen 
in many nearby early-types are driven by (low-level) merger-induced star formation 
within the last 3 Gyrs, contributing less than 10% of the stellar mass. If the ongoing 
mergers in this sample, which have a median mass ratio of 1:4, are representative of 
the nearby red merger population, then less than ^ 25% of the new stellar mass in the 
remnants is typically added through merger-induced star formation. While the dust 
extinction in the inter-stellar medium (ISM) in these galaxies is small {Eg^Jy < 0.1), 
the local dust content of the star-forming regions is, on average, a factor of ~ 3 higher. 
Finally, we use our theoretical machinery to provide a recipe for calculating the age of 
the most recent star formation event (^2) in nearby (z < 0.1) red early-type galaxies: 
log t2 [Gyrs] - 0.6±°°3 x [{NUV -u)-{g-z)- 1.73±°°3], where NUV, u, g and z are 
the observed photometric magnitudes of the galaxies in the GALEX/ SDSS filtersets. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters 



1 INTRODUCTION 

The formation histories of luminous, spheroidal galaxies 
have posed a serious challenge to models of galaxy forma- 
tion for the last few decades. The classical 'm onolithic col- 
lapse' hypothesis, that followed the model o f lEggen et al.l 
111962) for the formatio n of the Galaxy (e.g. lLarsorJll974l : 
IChiosi fc Carrardl2002l ). postulated that stellar populations 
in early-type galaxies form in short, highly efficient star- 
bursts at high redshift {z 3> 1) and evolve purely passively 
thereafter. Indeed, the optical properties of the early-type 
population and their strict obedience to simple scaling re- 
lations are remarkably consistent with such a simple for- 
mation scenario. The s mall scatter in the ear l y-type 'Fun- 
damental Plane' (e.g. Ijorgensen et al.l 119961 : ISaglia et al.l 
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1 19971) and its apparent lack of evo l ution with look-back time 
(e.g. iForbes et all Il998l: iPeebled I2OO2I : iFranxl Il993l . Il995l : 
Ivan Dokkum fc Franxlll996l ). the homogeneity and lack of 



redsh i ft evolution in their optical col o urs (e.g. Bower et al.l 
19921: iBendeiJ Il997l: IeuTs et all Il997l: IStanford et~all Il998l : 



Gladders et al.lll998l : lvan Dokkum et al.ll200Cl ) and evidence 
for short (< 1 Gyr) star formation timescales in these sys- 
tems, deduced from the over-abundance of a elements (e.g. 
iThomas et al.lll999l l. all suggest that the bulk of the stellar 
population in early-type galaxies does indeed form at high 
redshift (z > 2). 

Nevertheless, consensus, in recent years, has moved 
away from monolithic collapse towards a more gradual as- 
sembly of these systems in the framework of the stan- 
dard ACDM galaxy formation paradigm. Following the sem- 
inal work of iToomrd ([1977]), numerical simulations have 
consistently demonstrated that galaxy collisions end in 
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rapid merging - a fundamental feature of the standard 
model - and typically pr oduce spheroidal remnants (e.g. 
iBarnes fc Hernauistlll992l ) in cases where the mass ratios 
of the progenitors are large (> 1 : 3). Semi-analytical mod- 
els of galaxy formation, that incorporate these results, have 
had reasonable success in reproducing the morphological 
mix of the Universe and the (photometric) properties of the 



early - type population (e. g. Cole et a l. 2000; Hatton et al.l 



l2003l : iKavirai et al] l2005l : iDe Lucia et al. 200(^ 1. The pre 
dieted star formation histories (SFHs) of early-type galaxies 
(at least in clusters) in the semi-analytical framework have 
been shown to be quasi-monolithic, resulting in good agree- 
ment with the optical col ours and their evolution with red- 
shift (jKavirai et al.ll2005l ). As a result, the observed optical 
properties of the early-type population and their consistency 
with the expectations of monolithic collapse are not proof 
of a uniquely monolithic origin. 

A drawback of optical data is its relative insensitivity 
to small amounts of recent star formation (RSF). The op- 
tical spectrum remains largely unaffected by the minority 
of stellar mass that is expected to form in these systems 
at low and intermediate redshifts (z < 1), which makes it 
difficult to accurately measure early-type star formation his- 
tories (SFHs) over the last 8 billion years. However, while 
it does not impact the optical spectrum (within typical ob- 
servational and theoretical uncertainties) , a small mass frac- 
tion (a few percent) of young (a few tens of Myrs old) stars 
strongly affects the rest-frame UV spectrum shortward of 
3000A. Furthermore, the UV r emains largely unaffected by 
the age-metallicity degeneracy (jWorthevll 19941 ) that plagues 
optical analyses (Kaviraj et al. 2007a), making it an ideal 
photometric indicator of RSF. The advent of the GALEX 
space telescope (Martin et al. 2005) has provided us with an 
unprecedented opportunity to harness the sensitivity of the 
UV to young stars to quantify the presence of R SF in lo- 
cal ea rly-type galaxies. Following the early work of iYi et al.l 
(|2005l ). Kaviraj et al. (2007b; K07 hereafter), Schawinski et 
al. (2007a) and Kaviraj et al. (2008) have carefully stud- 
ied the UV propert i es of the early-type population at late 
epochs fsee lKavirail (|2008l ) for a review), to convert the UV 
flux observed in these systems into a measurement of the 
RSF, taking into account potential contributions to the UV 
spectrum from old, evolved stellar stages such as horizon- 
tal branch (HB) stars. It is worth noting that models where 
early-types host only old stars with realistic metallicity dis- 
tributions that satisfy both the mean metallicity and alpha- 
enhancements of nearby early-types clearly require RSF (see 
Section 6.2 in K07). These results indicate that early-types 
of all luminosities form stars over the lifetime of the Uni- 
verse, with massive systems accumulating 10-15% of their 
stars since z ^ 1, plausibl y through minor ra ergers with 
small, gas-rich companions (|Kavirai et al.ll2009l ). 

While such recent efforts have made a compelling case 
for RSF-driven UV excess in the early-type population, we 
should note some caveats to this analysis. The contamina- 
tion of the UV spectrum by old HB stars is a possibility 
at l ow redshift, e specially in giant elliptical galaxies (see 
e.g. iYi et al.lll997l . '1999). However, studies of the rest-frame 
UV at high redshift {z > 0.5), where the HB is absent, 
reveal a similarly large scatter in the rest-frame UV colour- 
magnitude relation (CMR) to what is observed in the nearby 
Universe, strongly suggesting that the principal drivers of 



the UV excess are not HB stars (Kaviraj et al. 2008). It is 
conceivable that there is some contribution to the UV flux 
from alternative hot stellar populations such as blue strag- 
glers - sta rs found above the turn-off in Hertzprung-Russell 
diagrams (|Wheelerl 19791 ) - typically within old globular clus- 
ters (GCs). Since they are like ly to form through stellar 
interactions IjMapelli et al.ll2006l ). these objects tend to in- 
habit the cores of GCs, where the stellar density is high. 
However, GCs (to which their constituent blue stragglers 
contribute a vanishing fraction of stellar mass) host a small 
fraction of the total stellar mass in the galaxy, making it 
unlikely that the blue straggler pop ulation are respons ible 
for the UV flux in elliptical galaxies (|Mochkovitchlll986l ) . 

We note that such alterative contributions to the UV 
flux in early-type galaxies have not been fully explored in the 
context of contemporary rest-frame UV data from GALEX 
at low redshift and deep optical surveys at high redshift. 
While RSF provides the most attractive solution, both in 
terms of the observed UV spectra of individual galaxies and 
in the context of the currently accepted galaxy formation 
paradigm, the (presumably minor) role of these alternate 
sources of UV flux remains an open question. A potentially 
robust discriminant between RSF and old stars is the spatial 
distribution of UV flux in the galaxy. If it is dominated by 
old stars (e.g. HB stars), the UV light is expected to be a 
smooth background, perhaps with a colour gradient (due to 
metallicity variations across the galaxy), but without large 
pixel-to-pixel variations or pronounced fine structure. A firm 
understanding of the origin of the UV flux would greatly 
benefit from high-resolution imaging of nearby early-types, 
e.g. using forthcoming instruments such as the HST-WFC3 
camera. 

In this work, we study the GALEX (UV) and SDSS (op- 
tical) properties of early-type galaxies on the optical red se- 
quence, originally studied bytvan Dokkum (2005, vDQ5 here- 
after). The novelty of the vD05 study was the unprecedented 
depth of the optical imaging employed to study the red 
galaxy population. Using ~27,000 second exposures on 4m 
class telescopes (equivalent to ~20 hrs on a 2.5m class tele- 
scope like the one used for the SDSS) vD05 found widespread 
signatures of recent merging in a large fraction of their early- 
type sample - more than 70% of bulge-dominated (E/SO) 
galaxies on the optical red sequence exhibit distorted mor- 
phologies, with predominantly red tidal features that extend 
to spatial scales > 50 kpc. The deep images allow a quan- 
titative measure of the perturbed morphologies of nearby 
early-type systems, which is not possible with the substan- 
tially shallower imaging offered by other surveys on 1-2. 5m 
telescopes. As an example, the effective exposure time of 
the SDSS - on which most of the recent early-type efforts 
are based - is only ~51s in each filter (York et al. 2000). 

GALEX UV imaging of the vD05 sample provides a 
unique opportunity to study the star formation history 
(SFH) of systems that are, to high confidence, recent merger 
remnants. We first establish that the vD05 sample is repre- 
sentative of the general luminous red galaxy population at 
similar redshifts. We then combine the GALEX/SDSS pho- 
tometry of the vD05 sample to quantify the RSF in these 
systems and explore the correlations between {NUV — r), 
the derived SFHs and the strength of the tidal distortions. 
Finally, we use our theoretical machinery to derive a simple 
recipe to translate the observed NUV and optical colours of 
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local early-type galaxies into an estimate of the age of the 
last star formation episode in these galaxies. 



2 THE GALAXY SAMPLE 

The original red galaxy sample of vD05 was selected 
from two deep extragalactic surveys: the Multi-wavelength 
Survey by Yale-Chile (MUSYC; Gawiser et al 2006) and the 
NOA O Deep Wide-Field Survey (NDWFS; Ijannuzi et all 
I1999I ). The combined dataset covers 10.5 deg^ of sky, reach- 
ing 1(7 surface brightness levels of /^ ~ 29 mag arcsec^^. 
We refer readers to Section 3 of vD05 for a comprehensive 
description of the dataset. 

For the purposes of this study the vD05 sample was 
cross-matched with both the SDSS DR4 and GALEX 
GR2q Out of the original 126 galaxies in the vD05 sam- 
ple, the cross-matching process yields 101 galaxies with both 
GALEX and SDSS photometry. The galaxies in this sample 
have a mean redshift of z ^ 0.11. The early-type galax- 
ies, that form the focus of this study, typically have super- 
L* luminosities (A/* ~ —21.15 in the nearby Universe; see 
Bernardi et al. 2003) and their masses range between lO^"* 
and lO"*^ Mq. 

Finally, we note that vD05 provides a quantitative mea- 
sure of the tidal distortions in each object by comparing 
galaxy images to elliptical galaxy models, generated using 
the ellipse task in IRAF (see Section 5.3 in vD05 for a de- 
scription of the procedure). This tidal parameter (denoted 
by T in this paper) measures the median absolute deviation 
of the (fractional) residuals from the model fit. Comparison 
of the value of T to visual inspection of the galaxy images 
indicates that, T < 0.08 in galaxies where the distortions 
are not obvious to the eye, 0.13 < T < 0.19 for weakly dis- 
turbed systems and T > 0.19 for systems that show strong 
tidal distortion. The maximum value of the tidal parameter 
in this sample is ~ 0.47. In the following sections, we use this 
tidal parameter to represent the extent of the morphological 
disturbance in each galaxy. 



3 EMISSION LINE PROPERTIES, SPATIAL 
DISTRIBUTIONS AND UV-OPTICAL 
COLOURS 

3.1 Emission line properties: AGN and 
star-formation activity 

The SDSS spectra of individual galaxies were used to 
determine the presence of Type II AGN from optical emis- 
sion line ratios, using a 'BPT'-type (e.g. Baldwin et al. 1981) 
analysis. 

In Figure [T] we present the BPT plot for the vD05 sam- 
ple. Grey dots indicate the AGN population in the low- 
redshift Universe, extracted by Kauffmann et al. (2003). 
Large circles denote galaxies in the vD05 sample. Galaxies 
are shown using filled circles if all four emission lines have 
S/N > 3. Galaxies where only [Nil] and Ha are detected 
are shown using open circles. Note that while the detection 
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Figure 1. A 'BPT'-type plot for the vD05 red galaxy sample, 
which indicates that only two galaxies lie near the star form- 
ing sequence and that the bulk of the emission-line activity is 
LINER-like. The grey dots indicate the AGN population in the 
low- redshift Universe, studied by Kauffmann et al. (2003). The 
large circles denote galaxies in the vD05 sample. Galaxies are 
shown using filled circles if all four lines have S/N > 3. Galaxies 
where only the [Nil] and Ha are detected are shown using open 
circles. 



of [Nil] and Hq allows us to classify these galaxies as AGN, 
the type of the AGN remains unclear, since the position on 
the y-axis depends on the strength of [QUI] and H/3 which 
are not robustly detected. 

It is remarkable that only two of the galaxies lie close 
to the 'star-forming' sequence of Kauffmann et al (which 
is demarcated by the solid curve in the BPT plot). This is 
strong evidence that the red mergers in vD05 are indeed 
mostly 'dry' and are not accompanied by s ignificant ongo- 
ing star formation (see also lWhitaker fc van Dokkurn 20081 ). 
The AGN fraction in this sample is ~ 11%, consistent within 
errors (~ ±3%) with the fraction found in the general early- 
type populat ion (15%) drawn from the entire SDSS at low 
redshift fe.g. lSchawinski et al.ll2007l ). It is also worth noting 
that most of the AGN in this sample are classified as LIN- 
ERs. A significant fraction of these systems could, in fact, be 
systems with declining star formation, where the 'LINER- 
like' emi^sionj£_djiyen_by post-AGB stars and white dwarfs 
fsee lStasihska et al.ll2008h . While further exploration of this 
issue is beyond the scope of this paper, we do not find com- 
pelling evidence of strong AGN activity in the vD05 sample. 



3.2 Spatial distributions and UV-optical colours 

We begin by studying some basic properties of the vD05 ob- 
jects. Figure [2] indicates the spatial coordinates (RA, DEC 
and redshift) of the galaxy sample (top) and its redshift 
distribution (bottom). Symbol sizes are proportional to the 
tidal parameter described in Section 2. Galaxies that con- 
tain active AGN are shown using open symbols and objects 
that are currently interacting are indicated using boxes. The 
redshifts of the galaxies in the top panel of this figure are 
shown colour-coded. We also indicate the angular sizes, at 
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z = 0.1 and z = 0.05, of a 3 Mpc structure (similar to the 
extent of the Coma cluster) . 

Given the high incidence of tidal features in the vD05 
early-type sample, we first check the spatial distribution of 
the objects to look for particular overdensities that might 
lead to the observed high frequency of interactions in the 
sample. The vD05 study did not consider the spatial distri- 
bution of the galaxy sample due to a lack of redshifts. They 
did estimate, however, that the median redshift of the galax- 
ies was z ~ 0.1, based on the optical {B — R) colour of an L* 
elliptical galaxy. Based on the spectroscopic redshifts from 
the SDSS (bottom panel of Figure [Sjl we find that the mean 
redshift of this sample to be ~ 0.11, in good agreement with 
the value estimated by vD05. 

We find that the sample occupies three main redshift 
peaks. Note, however, that the redshift binning used in Fig- 
ure[2]is 5z = 0.01 - which translates to a distance of ~50 Mpc 
at z = 0.1 in the standard cosmology - and that the peaks 
have widths of at least 5z = 0.02. Given the high accuracy 
of the SDSS redshifts (~ 10""*), the redshift distribution 
does not indicate that the galaxies are part of structures 
as compact as cluster cores (which are sub- Mpc structures). 
Furthermore, the transverse distances (in RA and DEC) are 
also quite large compared to, for example, the full extent of 
the Coma cluster (3 Mpc). Given that the centra l regions 
of ne arby clusters are dominated by red galaxies (JDresslei] 
[l980D, we do not find compelling evidence for the vD05 sam- 
ple to be drawn out of a small number of overdense regions. 
The galaxies therefore appear to be, on average, in the 'field'. 

The top panel of Figure |3] shows the UV and optical 
colour-magnitude relations for the vD05 objects. Small black 
dots indicate the nearby (0 < 2; < 0.1), massive (> L*) 
early-type population, drawn from the SDSS DR4 following 
the method used by K07. Briefly, early-type galaxies are se- 
lected from the SDSS using the pipeline parameter f racdev, 
which provides a probabilistic measure of how well a 'de Vau- 
couleurs' proflle fits the galaxy's light profile. Galaxies with 
f racdev > 0.95 are extracted and then visually inspected to 
remove contaminants such as face-on spirals. Galaxies blue- 
ward of [NUV — r) ~ 5.5 (indicated by the green horizontal 
line) are very likely to have had some star formation within 
the last Gyr (see Section 3 in K07). It is clear that, while 
the galaxies in this sample are optically red, they exhibit a 
wide spread in UV colours. 

We note that a central AGN could contaminate the UV 
spectrum of a galaxy, perhaps enhancing the blue UV colour. 
However, the contamination from a Type II AGN is likely 
to be less than around 15% in UV flux, which translates to 
around 0.15 mags in the [NUV — r) colour, much smaller 
than the spread in the UV CMR (~ 6 mags; see Salim et al. 
2007). Typically, blue early-type colours {NUV - r < 5.5) 
are not restricted to galaxies hosting AGN, either in the 
sample studied here or in the general early-type population 
studied by K07. Furthermore, the quality of the SED fitting 
that is used to calculate star formation histories (see the 
next section) is equally good in galaxies which carry BPT 
signatures of AGN and those that do not show any signs 
of AGN, indicating that there is no measurable contribu- 
tion from a power-law component. Finally, a study of the 
GALEX images of nearby AGN hosts indicates that the UV 
emission is extended, making it unlikely that it comes from 
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Figure 2. TOP: Spatial coordinates (RA, DEC and redshift) of 
the galaxy sample studied here. BOTTOM: Comparison of the 
redshift distribution of the general galaxy sample to that of only 
the interacting population. Our results indicate that the vD05 
galaxies are generally in the 'field'. Symbol sizes arc proportional 
to the tidality parameter described in Section 2. Galaxies with 
active AGN are shown using open symbols and objects that are 
currently interacting arc indicated using boxes. Open boxes in- 
dicate ongoing mergers with AGN activity. We also indicate the 
angular sizes, at 2 = 0.1 and z = 0.05, of a 3 Mpc structure (sim- 
ilar to the extent of the Coma cluster). The redshift binning in 
this plot is 5z ~ 0.01, which translates to a distance of ~50 Mpc 
at 2 = 0.1 in the standard cosmology. 



a central source (Kauffmann et al. 2007). Thus, the presence 
of Type II AGN does not affect the analysis in this paper. 

We begin by comparing the vD05 galaxies to the general 
early-type population. For a direct comparison we restrict 
the analysis to the K07 optical red sequence {g — r > 0.65) 
in the redshift range 0.05 < 2: < 0.1 and compare the UV 
colour distribution of the two samples (bottom panel of Fig- 
ure[3|. We find that the UV distributions are consistent with 
each other (within counting errors), although the vD05 sam- 
ple does show a slight excess of objects in two intermediate 
colour bins around {NUV - r) ~ 4.3 and {NUV - r) ~ 4.7. 

Recall, however, that the method for selecting the gen- 
eral early-type galaxy population in K07 involves extracting 
galaxies whose light profiles in the g, r and i-bands match 
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Figure 3. TOP: The optical and UV colour-magnitude relations 
(CMRs). While selected to be optically red, the vD05 sample has a 
large range of UV colours. Red symbols indicate early-type galax- 
ies while blue symbols indicate spirals. Symbol sizes are propor- 
tional to the tidal parameter described in Section 2. Galaxies that 
contain active AGN are shown using open symbols and objects 
that are currently interacting are indicated using boxes. Small 
black dots indicate the nearby (0 < z < 0.1), massive (> L*) 
early-type population drawn from the SDSS following the method 
used by K07. Galaxies blueward of {NUV — r) ~ 5.5 (indicated 
using the green horizontal line) are very likely to have had some 
star formation within the last Gyr (see Section 3 in K07). BOT- 
TOM: Comparison between the UV colours of the general early- 
type population on the optical red sequence {g — r > 0.65) and 
the vD05 sample. See text for more discussion. 



a 'de Vaucouleurs' profile very closely (by setting f racdev> 
0.95 when selecting objects from the SDSS database). To 
quantify the effect of this criterion on the colour compar- 
ison, we study the variation of the {NUV — r) colour in 
the vD05 sample with their fracdev values (inset in the 
bottom panel of Figure [3]). About 57% of the vD05 galax- 
ies with fracdev< 0.95 have {NUV - r) < 5. The frac- 
tion of early-types that were missed from K07 as a result of 
the f racdev> 0.95 criterion is around 20%, from early-type 
catalogs selected by eye-inspection alone (i.e. no fracdev 
cut; Kevin Schawinski priv. comm.). Assuming a red se- 
quence fraction of ^ 70% (since K07 estimated that 30% 
of the early-type population were blue) implies an ~ 8% 
(57% X 20% X 70%) increase in the fraction of general early- 
types with colours blueward of NUV — r ~ 5. 

This correction removes most of the discrepancy in the 
UV colour distributions. However, the fraction of objects in 
the colour bin around {NUV — r) ~ 4.3 remains inconsistent 
within one-sigma uncertainties, although the histograms are 
consistent at the two-sigma level. A s imple Kolmogorov- 
Smirnov (KS) test (see e.g. I Walll 119961 ) indicates that the 
probability that the two populations are drawn from the 
same parent distribution is ~72%. A similar conclusion is 
drawn if we use a Kuiper test (which maintains sensitiv- 
ity near the tails). Thus, we find good overall agreement 
between the colour distributions, with the difference plau- 
sibly caused by cosmic variance, since the sky coverage of 
the SDSS DR4 (from which the early-type population is de- 
rived) is more than two orders of magnitude larger than that 
of the vD05 sample. 



3.3 Correlation between colour and tidal 
distortion 

In Figure |4] we plot the UV and optical colours against the 
tidal parameter. Note that the vertical scale in the top panel 
of this figure, which shows the {NUV — r) colour, spans 4 
mags, while that for the {g — r) colour spans only 0.3 mag. 
We explore the strength of the correlation between colour 
and tidal distortion using Spearman's rank correlation test, 
which produces a non-parametric measure (p) of the statisti- 
cal dependence between two variables. The value of p varies 
between -1-1 and -1, which indicate perfect positive and neg- 
ative correlations respectively, p ~ indicates a complete 
lack of correlation between the variables. 

We find that, for the relationship between {NUV — r) 
and the tidal parameter, p — —0.43. The significance of 
this value is 0.0032. This is the probability of finding a 
correlation between {NUV — r) and the tidal parameter 
purely through statistical accident. Since this number is very 
small, this implies that the correlation is statistically sig- 
nificant. The corresponding values for the {g — r) colour 
are p — —0.08, with a significance of 0.35. Note that the 
probability of a chance correlation is higher for the optical 
colour essentially due to the correlation being weaker. We 
find, therefore, that a measurable correlation exists between 
both the optical and UV colours and tidal distortion in the 
galaxy. Not unexpectedly, the correlation is much stronger 
(and has higher confidence) in the UV colour, since it is 
much more sensitive to small amounts of RSF. 

The reasonably large scatter in the relationship between 
colour and tidal distortion can be understood by noting that 
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Figure 4. Tidal parameter plotted against {NUV — r) and (g — r) 
colours. A statistically significant correlation exists between the 
tidal parameter and both the UV and optical colours (see text in 
Section 3.3 for details). The symbol size scales with the size of the 
tidal distortion. Galaxies that host Type II AGN are shown using 
open symbols. Bluer galaxies show higher tidal distortions, with 
the trend being much stronger in the {NUV — r) colour, due to 
its higher sensitivity to small amounts of recent star formation. 



the stellar populations that are responsible for the tidal fea- 
tures and those that drive the UV colours tend to be dif- 
ferent. Assuming that the RSF is driven by merging (see 
also Kaviraj et al. 2010), numerical simulations (e.g. Kavi- 
raj et al. 2009) indicate that the tidal features form from 
the underlying stellar populations of the accreted satellite, 
while the UV colour is sensitive to the gas content of the 
satellite and how far star formation has progressed when 
the remnant is viewed. For example, two satellites with sim- 
ilar underlying stellar populations but diflerent gas contents 
might yield similar tidal distortions (if viewed at the same 
point during the merger process) but exhibit different UV 
colours. Thus, while a high tidal distortion should, in gen- 
eral, be expected to coincide with a blue UV colour because 
the remnant is viewed at an earlier phase in the merger/star 
formation process and the stellar populations are younger, 
the spread in the satellite gas contents will induce a scatter 
into the correlation. 



4 RECENT STAR FORMATION HISTORIES 

4.1 Parameter estimation 

A key objective of this work is to quantify the RSF in the 
vD05 sample. We are primarily interested in exploring the 
age of the last star formation event (which is presumably 
linked to the morphological disturbances observed in the 
deep optical images) and the stellar mass fractions con- 
tributed by these events. 

We estimate parameters governing the star forma- 
tion history (SFH) of each galaxy by comparing its 
(FUV, NUV, u, g, r, i, z) photometry to a library of synthetic 
photometry, generated using a large collection of model 
SFHs, specifically optimized for studying spheroidal galax- 
ies at low redshift. As we describe below, our scheme decou- 
ples the RSF episode from the star formation that creates 
the bulk, underlying population. We choose a parametrisa- 
tion for the model SFHs that both minimises the number 
of free parameters and captures the macroscopic elements 
of the star formation history of spheroidal galaxies in the 
low-redshift Universe. 

Since the underlying stellar mass in spheroidals forms 
at high redshift and over short timescales, we model the 
underlying stellar population using an instantaneous burst 
at high redshift. We put this first (primary) instantaneous 
burst at z=3. Note that changing this to z=2, or even z=l, 
does not affect our conclusions about the RSF, because the 
first burst does not contribute to the UV. A large body of 
recent evidence suggests that the star formation in these 
syst ems in the nearby U niverse is driven by minor mergers 
(see lKavirai et al.ll201(J . and references therein). This star 
formation is bursty and we model the RSF episode using 
a second instantaneous burst, which is allowed to vary in 
age between 0.001 Gyrs and the look-back time correspond- 
ing to z = 3 in the rest-frame of the galaxy, and in mass 
fraction between and 1. Our parametrisation is similar to 
previous o nes used to study elli ptical galaxies at low red- 
shifts (e.g. iFerreras fc Silkll2000l ). We briefly note that our 
scheme differs in construction from other recent work that 
explore the general galaxy population at higher redshift s, 
without reference to morphology (e.g. iMuzzin et al.ll2009l ). 
In particular, the decoupling employed in this study between 
RSF and the bulk stellar population is difficult to apply at 
high redshift, where the underlying stellar populations are 
actively being assembled. Thus, our particular parametrisa- 
tion works best in spheroidal galaxies at low redshift and is 
not directly comparable to more general schemes, especially 
in the high redshift regime. 

To build the library of synthetic photometry, each 
model SFH is combined with a single metallicity in the range 
O.IZq to 2.5Z0 and a value of dust extinction parametrised 
by Eb-v in the range to 0.5. The dust model empl oyed in 
this st udy is the empirical dust prescription of Calzc tti et al.l 
(|200d ). Photometric predictions are generated by combin- 
ing each model SFH with the chosen metallicity and Eb-v 
values and convolving with the stellar models of lYJ (|2003l ) 
through the GALEX and SDSS flitersets. The model library 
contains ~ 750, 000 individual models. Note that the FUV 
and NUV filters have effective wavelengths of ^ 1500A and 
~ 2300A respectively. 

Since our galaxy sample spans a range of redshifts, 
equivalent libraries are constructed at redshift intervals of 
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Sz — 0.01. A fine redshift grid is essential in such a low 
redshift study because a small change in redshift produces 
a relatively large change in look-back time over which the 
UV flux can change substantially, inducing 'K-correction- 
like' errors into the analysis. 

The primary free parameters in this analysis are the age 
(f2) and mass fraction (/2) of the second burst (the mass 
fraction of the primary burst is simply 1 — /2). Secondary 
parameters of interest are the dust properties of the system. 
In particular, the dust content of the primary burst, Eg^y, 
(which is effectively tracing the 'underlying' population of 
the galaxy) should provide a measure of the extinction in the 
inter-stellar medium (ISM), while the dust in the secondary 
burst (Eg^y) should reflect the dust in the star formation 
regions and their vicinity. 

In each case, the value of the free parameters are esti- 
mated by comparing each observed galaxy to every model 
in the synthetic library, with the likelihood of each model 
(exp— x'^/S) calculated using the value of x^j computed in 
the standard way. In addition to the observational errors 
for each object, we assume additional uncertainties that ac- 
count for uncertainties in the models (e.g. Yi 2003), model 
offsets from observational data (e.g. Eisenstein et al. 2001; 
Maraston et al. 2009) and systematic errors in GALEX and 
SDSS photometry (see e.g. Ivezic et al. 2004; Morrissey et 
al. 2007). These are taken to be 0.05 mags for the optical 
filters, 0.2 mags in FUV and 0.1 mags for the NUV pass- 
band. The two types of uncertainty are added in quadrature. 
From the joint probability distribution, each parameter is 
marginalised to extract its one-dimensional probability den- 
sity function (PDF). We take the median of this PDF as 
the best estimate of the parameter in question and the 16 
and 84 percentile values as the 'one-sigma' uncertainties on 
this estimate. In the analysis that follows we present these 
median parameter values. 

The leverage in f2 and the quality of the t2 fits depends 
critically on our access to the rest-frame UV, which hosts 
most of the flux from hot, youn g main sequence stars. Note 
that IWhitaker fc van DokkumI ()2008l l have demonstrated 
that the spheroidal galaxies in the vD05 sample are vir- 
tually dustless, so that there should be very little repro- 
cessing of the UV-o ptical light into the infrared passbands. 
iMuzzin et al.l l|2009h have shown that the inclusion of rest- 
frame near-infrared (NIR) filters does not alter the mean 
values of fitted parameters (although they do reduce the de- 
rived uncertainties) for galaxies in the high-redshift {z ~ 2) 
Universe, which should be significantly dustier than our low- 
redshift spheroidal sample. Since our sample is dust-poor, it 
seems reasonable, therefore, to suggest that the addition of 
NIR filters will leave our conclusions regarding the derived 
parameters unchanged. Finally, a similar analysis on a small 
sample of SDSS early-type galaxies which have, in addition 
to GALEX, NUV photometry from the XMM 'Optical Mon- 
itor (OM)' (the OM has three NUV filters around GALEX 
NUV, see Mason et al. 2001) indicates that a larger num- 
ber of UV filters leaves the median values of i2 and /2 un- 
changed, with the uncertainties decreasing by ~10%. Thus 
the size of our filterset, both in the UV and in the optical 
spectrum, are optimal for the aims of this paper. 



4.2 Star formation history parameters 

4-2.1 Age and mass fraction of the second burst 

Figure [5] shows the age (ti) and mass fraction (/2) of the 
RSF estimated to have formed in the second burst in each 
of the vD05 early-types. Galaxies are colour-coded based on 
their {NUV — r) colour and average uncertainties in the 
parameters are indicated in the plot. Note that, although 
the U V colour changes rapidly with the value of t2 , the same 
colour can be degenerate with a wide range of mass fractions 
/2 (see Figure 1 in Kaviraj 2008). Since this effect becomes 
more pronounced as the starburst ages, the mass fraction 
errors become larger at higher ages. As a result, while we 
can typically constrain the age of the second burst with 
good accuracy, the uncertainty in the mass fraction could 
span almost a decade. 

We find that the bulk of the objects in our sample have 
experienced a star formation event in the last 3 Gyrs. Con- 
sidering galaxies which have well-constrained mass fraction 
uncertainties, the median values of /2 are typically less than 
10%. 5% of the early- type sample has experienced a star- 
formation event within the last 0.1 Gyrs. The very low mass 
fractions in this subset may indicate that the star formation 
in these events is still ongoing. 

Note that the median mass ratio of the ongoing mergers 
in the vD05 sample is 1:4. If the ongoing mergers are repre- 
sentative of the red merger population, then less than 25% 
of the new stellar mass in the remnants is contributed by 
star formation. The rest is 'dry accretion' of existing stars 
in the progenitor galaxies. Note, however, that the errors in 
the mass fraction estimates are typically quite large and this 
statement probably carries with it some uncertainty. Using 
the lower limit of the error bars in the mass fraction would 
suggest a lower fraction of less than 10% being contributed 
by star formation, which is consistent with the recent liter- 
ature fe.g. Ivan Dokkum et al.ll2010r i. 

In Figure [6l we plot the age (fe) and mass fraction (/2) 
of the second burst against the tidal parameter. Recall that 
t2 is effectively the look-back time to the last star formation 
event in the galaxy. We find that galaxies with larger tidal 
distortions show lower values of i2, consistent with the ex- 
pectation that the intensity of the tidal distortion decreases 
with time, resulting in the observed correlation between the 
tidal parameter and t2. The best-fit relation is described by: 



logT^ 



-o.gr 



Q2i±005 xlogt2 [Gyrs] 



(1) 



A weaker correlation is found between the mass fraction 
formed in the recent starburst (/2) and the tidal parameter, 
as shown in the bottom panel of Figure |6l Note that the 
range on the x-axis in this plot is much larger than that in 
the t2 vs tidal parameter plot (top panel). Since the tidal 
features are composed of the the underlying stellar mass 
of the accreted satellite and not the young stars emerging 
from the starburst (which tend to form preferentially in the 
central regions of the remnant), the tidal distortion should 
be relatively insensitive to /2. However, a residual correla- 
tion should be expected, since a larger tidal distortion indi- 
cates that the merger is being observed at an earlier stage 
in the relaxation process. Hence the mass fraction formed is 
smaller than what it would be if the same remnant was ob- 
served later, when the tidal distortion would have decreased 
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Figure 5. Ages (i2) and mass fractions (/2) of the second star 
formation episode in the early-type galaxies studied here. The 
bulk of the vD05 sample has experienced a star formation event 
in the last 3 Gyrs, which typically adds less than 10% to the 
stellar mass of the galaxy. Symbol sizes are proportional to the 
tidal parameter described in Section 2. Galaxies that contain ac- 
tive AGN are shown using open symbols and objects that are 
currently interacting are indicated using boxes. The galaxies are 
shown colour-coded based on their (NUV — r) colour. 



and a larger portion of the star formation would have been 
completed. 

The observed trends between colour, i2, /2 and tidal 
distortion can be explained by the fact that galaxies with 
higher tidal distortions are at an earlier stage in the merger 
process and the star formation episode that is induced by 
it. As time passes the remnant relaxes and the star forma- 
tion activity induced by the merger decreases, as the gas 
fuelling this process gradually runs out. This results both in 
lower tidal distortions, redder colours, and increasing mass 
fractions, accounting for the trends observed in this study. 

We note that our results are in good agreement with the 
recent study of Sanchez-Blazquez et al. (2009), who studied 
the vD05 sample using spectral line diagnostics. 




1 



Figure 6. The age (42) and mass fraction (/2) of the recent star 
formation in the vD05 sample, plotted against the strength of 
their tidal distortions. Galaxies that have commenced their star 
formation more recently (i.e. lower values of ^2) show larger tidal 
distortions. A weak trend also exists between the mass fraction /2 
and the tidal parameter, in the sense that smaller mass fractions 
correspond to higher values of the tidal parameter. 



4-2.2 Metallicity and dust properties 

The median metallicity values for the galaxies in our sample 
are in the range 0.8 to 1.4 Zq, as might be expected for 
spheroidal galaxies. Recall, from Section 4.1, that a single 
metallicity is fitted to the entire galaxy, so that the derived 
metallicity is a represenative value for the combined stellar 
population (old + young) in the galaxy. 

In Figure [7] we present the dust properties of the galax- 
ies in our sample. We find that the ISM extinction in these 
objects, which is traced by the dust content of the pri- 
mary burst (Eg^y), is generally less than 0.1, consistent 
with the expected properties of early-type galaxies as rel- 
at ively dust-poor systems. In agr eement with the findings 
o n Whitaker & van Dokkum (20081) . we find that the redness 
of the vD05 objects is not due to large amounts of dust but 
simply because the bulk of their stellar populations is old. 

While the macroscopic ISM seems to be relatively dust- 
free, we typically derive larger values for Eg^y (which is rep- 



resentative of the dust in star-forming regions in the galaxy) . 
The bulk of the sample satisfies Eg'iy < 0.4. Figure [7] in- 
dicates that the values of Eg'iy can be a few factors larger 
than Eg^y, which is consistent with the fact that the extinc- 
tion around star forming regions (specifically in molecular 
clouds) is expected t o be several times lar ger than that due 
to the ISM alone (e.g. lCharlot fc FallllJOOol ) . The mean value 
of Eg^y/Eg'iy Is ~ 3 (uote that the theoretical estima te, 
for starburst galaxies, derived by ICharlot fc Falll (|200G| ) is 
also ~ 3). While E%'^y applies to only a minority (< 10%) 
of the stellar mass in the galaxy, it is worth noting that the 
higher extinction appears to affect the second burst even 
when t2 is substantially greater than the age s of molecu- 
lar clouds , which are typica lly 10-30 Myrs (e.g lBlitz fc Shul 
|l980; Har tmann et al.ll200ll ). This may indicate that regions 
of the galaxy where the RSF originates remain dusty. In ad- 
dition, the parametrisation used in this study does not take 
into account the timescale over which the RSF takes place 
(since all bursts are assumed to be instantaneous). In real- 
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Figure 7. The dust content of the primary burst plotted again 
the dust content of the secondary burst. The dust extinction of 
the primary burst (which represents the ISM) is typically less 
than 0.1, consistent with early-types being dust-poor. The extinc- 
tion in star-forming regions is, on average, a factor of 3 higher. 
Symbol sizes arc proportional to the tidality parameter described 
in Section 2. Galaxies that contain active AGN arc shown using 
open symbols and objects that arc currently interacting are indi- 
cated using boxes. The galaxies are colour-coded based on their 
{NUV - r) colour. 



ity, only the youngest stars in the RSF age distribution will 
probably suffer high extinction but the lack of a timescale in 
our analysis leads to the entire second burst suffering higher 
extinction to account for the plausible effect of molecular 
clouds. 



4.3 A prescription to estimate the age of the last 
star formation event 

The parameter estimation performed in this study allows us 
to develop a recipe to estimate the age of the most recent star 
formation event, given the UV and optical photometry of the 
galaxy in question. We focus on a recipe of the age rather 
than the mass fraction because the smaller uncertainties in 
the age make the age estimation significantly more robust. 

We note first that the median value for i2 used in 
this analysis is extracted from its marginalised PDF, so 
that the effects of all other parameters (mass fraction, dust 
and metallicity) have been integrated out. In the top panel 
of Figure |8l we plot the median value of t2 against the 
{NUV — r) colour. We also show the equivalent plot for 
the optical {g — r) colour in the inset. 

While a clear correlation exists between i2 and the 
{NUV — r) colour, there is no equivalent trend with the 
optical colour. This plot demonstrates both the usefulness 
of the UV spectrum in quantifying the age of the RSF in 
early-type galaxies and the insensitivity of the optical spec- 
trum to the residual star formation in these objects. 

Although a correlation exists between i2 and {NUV—r), 
we note that the {NUV — r) colo ur is dust sensitive . 
Standard dus t exti nction laws (e.g. ICalzetti et al.l I200C1 : 
ICardelli et al.l 119891 ) indicate that the extinction in the 
{NUV ~ r) colour, A^^^"'' ~ 6 x Eg-v- Recalling that 
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Figure 8. TOP: The {NUV — r) colour plotted against the age 
of the second burst. BOTTOM: The relatively 'dust-insensitive' 
{NUV — u) — {g — z) double colour plotted against the age of the 
second burst. The solid line indicates a linear least-squares fit to 
the relation, which allows us to estimate the age of the last star 
formation event in early-type galaxies directly from their observed 
NUV, u, g and 2-band photometry. 



the estimate for t2 is marginalised over all other parameters 
(including dust), we would, ideally, like to compare it to an 
observational quantity which is also insensitive to dust. 

Since the extinctions in colours are linear combinations 
of the extinctions in the constituent filters, we construct a 
double colour, {NUV — u) — {g — z), that minimises the 
dust sensitivity, given the filter combinations available to 
us for the bulk of t he vD05 samp l e. Th e extinction in this 
colour (assuming a lCalzetti et al.l l|200Cl ) extinction law) is 
^(Nuv-u)-{g-z) ^ -^5 ^ Eb-v, a factor of 4 less than that 
in {NUV — r). Not unexpectedly, plotting t2 against {NUV ~ 
u) — {g — z) produces the same quantitative trend but with 
less scatter (bottom panel of Figure |8]). The best- fit line 
(shown using the solid red line) yields: 



logi2 [Gyr 



0.6^ 



X [{NUV -u)-{g-z)-1.7r 



(2) 

Eqn (2) allows us to estimate the age of the last star 
formation event in early-type galaxies, directly from their 
observed (i.e. not K-corrected) NUV, u, g and z-band pho- 
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tometry. We note that there appears to be a possible flat- 
tening of the slope in the correlation for the very bluest 
{NUV - r < 3.8 or {NUV -u)-{g-z) <0) objects (which 
have been excluded from the fitting above). These objects 
may obey a different scaling, although the small number of 
galaxies in this part of the parameter space makes it difficult 
to derive a robust conclusion. It is possible that, given the 
very small values of i2 , the star formation activity in these 
events is in its very early stages, displacing them from the 
main locus. 



5 CONCLUSIONS 

We have combined GALEX (UV) and SDSS (optical) pho- 
tometry to derive the recent star formation histories of a 
sample of ~ 100 red sequence galaxies, many of which ex- 
hibit widespread signs of disturbed morphologies (in deep 
optical imaging) that are consistent with recent merging 
events. Originally studied by van Dokkum (2005; vD05 here- 
after), more than 70% of bulge-dominated galaxies in this 
sample show tidal features at a surface brightness limit of 
/i ~ 28 mag arcsec"^. An analysis of the spatial (RA, DEC 
and redshift) distributions of the vD05 objects indicates that 
the sample is not drawn from overdense regions such as 
cluster cores and can be considered to be part of the gen- 
eral 'field'. While the vD05 galaxies are optically red, they 
show a wide dispersion in UV colours that is attributable 
to low level RSF, akin to what has been found recently in 
the general early-type population. Comparison of the vD05 
UV colour distribution to that of the general early-type red 
sequence (drawn from the SDSS) using a KS test indicates 
that they are drawn from the same parent distribution. 

A statistically significant correlation is found between 
galax;y colour and the intensity of the tidal distortion, in 
the sense that redder objects show smaller tidal distortions. 
A stronger trend is observed in the UV colour, since it is 
more sensitive to low-level recent star formation. Not unex- 
pectedly, a correlation also exists between the parametrised 
look-back time to the recent star formation event (^2) and 
the tidal distortion, in the sense that the tidal distortion de- 
creases with increasing values of f2. Only a weak trend exists 
between the mass fraction formed in the event (/2) and the 
tidal parameter, with larger values of the tidal parameter 
corresponding to smaller values of /2. The observed trends 
can be explained by the fact that galaxies with higher tidal 
distortions are at an earlier stage in the merger process and 
the star formation episode that is induced by it. As time 
passes the remnant relaxes and the star formation activity 
decreases, as the gas fuelling this process gradually runs out. 
This results in lower tidal distortions, redder colours and in- 
creasing mass fractions, as found in the trends observed in 
this study. 

The bulk of the vD05 sample has experienced a star for- 
mation event in the last 3 Gyrs. The median values of the 
mass fractions are typically less than 10%. 5% of the sample 
has experienced a star-formation event within the last 0.1 
Gyrs and the star formation episodes in these objects ap- 
pear to be ongoing since the mass fractions formed in these 
events is very small (< 0.1%). If the ongoing mergers in the 
vD05 sample, that appear to have a median mass ratio of 1:4 
(see Section 6.2 in vD05), are representative of nearby red 



mergers, then less than ~25% of the new stellar mass in the 
remnants is contributed by star formation. The rest is 'dry 
accretion' of existing stars. We should note, however, that 
the large uncertainties on the mass fraction estimates (see 
Figure [5]) makes this statement somewhat uncertain and us- 
ing the lower limit of the error bars on the /2 values suggests 
a lower fraction of less than 10% being contributed by star 
formation. 

The parameter estimation performed in this study al- 
lows the primary burst (which takes place at z — 3 and 
contributes the underlying stellar population) and the sec- 
ondary burst (which represents the most recent burst of star 
formation and dominates the UV flux) to have different dust 
extinction values. We find that the dust extinctions in the 
primary burst {Eg^y) are typically less than 0.1, consis- 
tent with early-type galaxies having ISMs that are dust- 
poor. However, while the macroscopic ISM seems to be rel- 
atively dust-free, we typically derive larger values for the 
dust extinction in the second burst (Eg^y), with the bulk 
of the sample satisfying Eg^y < 0.4. The mean value of 
Eb'^v/Eb-v is ~ 3. We note, however, that, while Eg'ty 
applies to only a minority (< 10%) of the stellar mass in 
the galaxy, the higher extinction appears to affect the sec- 
ond burst even when its age (^2) is substantially greater 
than the ages of molecular clouds, which are typically 10-30 
Myrs. This may indicate that regions of the galaxy where 
the RSF originates remains dusty. It could also be a result 
of the parametrisation used here, which does not take into 
account the timescale over which the RSF takes place (since 
all bursts are assumed to be instantaneous). In reality, only 
the youngest stars in the RSF age distribution will prob- 
ably suffer high extinction but the lack of a timescale in 
our analysis leads to the entire second burst suffering higher 
extinction to account for the plausible effect of molecular 
clouds. 

We have used our analysis to develop a recipe to esti- 
mate the age of the most recent star formation event, given 
the observed UV and optical photometry of the galaxy in 
question. We find a well-defined correlation between the age 
of the second burst (^2), which is marginalised over all other 
free parameters and the relatively dust-insensitive double 
colour, (NUV ~ u) — (g — z). Eqn (2) allows us to esti- 
mate the look-back time to the last star formation event 
in red early-type galaxies, directly from their observed (i.e. 
not K-corrected) NUV, u, g and z-hand photometry. The 
slope in the correlation appears to fiatten for the very bluest 
{NUV~r < 3.8 or (NUV-u)-{g-z) < 0) objects, possibly 
due to the RSF being ongoing in these systems. 

Since the vD05 sample appears to be representative of 
the general field population of spheroidal galaxies in the 
nearby Universe (Section 3.2), the results of this study sug- 
gest that the widespread low-level star formation that has 
been reported in nearby spheroidals is driven by recent 
mergers. The vast majority of these events do not involve 
large amounts of gas, indicating that they are either dry ma- 
jor mergers or minor mergers, where the spheroid accretes a 
small gas-rich satellite (see also Kaviraj et al. 2010). Never- 
theless, the spheroidal galaxy population, at present-day, is 
constantly evolving, both dynamically through interaction 
with companions and in terms of low-level star formation 
activity induced by these interactions. 
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